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The presence of spin-orbit scattering within an aluminum nanoparticle affects measurements 
of the discrete energy levels within the particle by (1) reducing the effective g-factor below the 
free-electron value of 2, (2) causing avoided crossings as a function of magnetic field between 
predominantly-spin-up and predominantly-spin-down levels, and (3) introducing magnetic-field- 
dependent changes in the amount of current transported by the tunneling resonances. All three 
effects can be understood in a unified fashion by considering a simple Hamiltonian. Spin-orbit scat- 
tering from 4% gold impurities in superconducting aluminum nanoparticles produces no dramatic 
effect on the superconducting gap at zero magnetic field, but we argue that it does modify the nature 
of the superconducting transition in a magnetic field. 
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I. Introduction 

For decades, systematic studies of the quantum- 
mechanical energy levels of atoms and atomic nuclei have 
provided an understanding of the forces governing these 
systems. Recently, it has also become possible to mea- 
sure the discrete "electrons-in-a-box" energy levels within 
semiconductor quantum dots and metal nanoparticles 
[^|-^. Experiments have shown that different classes 
of forces and interactions acting on the electrons inside 
these materials affect the level spectra in distinguishable 
ways. Therefore, just as in atomic and nuclear physics, 
the discrete spectra in these condensed matter systems 
can provide a tool for understanding the interactions 
which influence electronic structure, uncovering effects 
that are not clearly visible if the individual quantum lev- 
els in the system cannot be resolved. The consequences of 
superconducting pairing interactions 1 4 and more gen- 
eral electron-electron interactions |^,7| have previously 
been analyzed for the case of aluminum nanoparticles. 
In this report, we discuss spin-orbit (S-0) interactions, 
resulting both from accidental defects in the Al nanopar- 
ticles and from gold dopants. We examine how S-0 scat- 
tering affects both the energies of the quantum levels and 
the amount of tunnel current which may be carried by 
each state. We find that the magnetic-field dependence 
of these quantities may be understood in a unified fashion 
within a simple model. The effects of S-0 scattering on 
the superconducting properties of an aluminum nanopar- 
ticle are also discussed. 

The study of S-0 scattering within metals has a long 
history. The metal samples of the types traditionally 



examined are large enough that the electronic states ef- 
fectively form a continuum. In this case, the quantity of 
primary experimental interest in S-0 studies is the rate 
at which S-0 interactions cause the spin of an electron as- 
sumed to be initially in a pure spin- up or spin-down state 
to be scattered into continuum states with opposite spin. 
This rate can be measured using weak localization exper- 
iments for disordered metal samples or, alternatively, 
by analyzing the form of the spin-dependent density of 
states determined by tunneling between thin supercon- 
ducting films in a parallel magnetic field S-0 interac- 
tions are of fundamental theoretical importance because 
their presence changes the symmetry properties of the 
Hamiltonian. For instance, the statistics of the energy 
levels in chaotic time-reversal-symmetric quantum dots 
are predicted to change from the orthogonal distribution 
in the absence of S-0 scattering to the symplectic distri- 
bution for a strong S-0 interaction, with a corresponding 
increase in the strength of the effective energy-level re- 
pulsion [ p^|JTl| ] . Perhaps the most dramatic consequences 
of S-0 coupling in metals occur in ferromagnets, since 
the S-0 interaction underlies the phenomena of magnetic 
anisotropy and the anomalous Hall effect. 

An analysis of the effects of S-0 interactions in metal 
nanoparticles requires a somewhat different viewpoint 
than for larger devices with a continuum density of states. 
Considering basic symmetries, the Hamiltonian operator 
describing electrons within a metal sample does not com- 
mute with the components of the total electronic spin op- 
erator in the presence of the S-0 interaction. This means 
that it is not possible to construct a set of basis states 
which are simultaneously eigenstates of both the energy 
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and Sz- The discrete energy eigenstates, through which 
electron tunneUng occurs in a metal nanoparticle, will 
thus necessarily be linear superpositions of pure spin-up 
and pure spin-down states, with the extent of admixture 
determined by the magnitude of S-0 matrix elements. 
Because these discrete energy eigenstates defined in the 
presence of the S-0 interaction are in fact well-defined 
energy eigenstates, the S-0 interaction does not lead to 
any decrease in their lifetime. For this reason, the ex- 
perimental quantities of interest in this paper will not 
be scattering rates, but rather shifts in the energies of 
the electronic states and changes in the tunneling cur- 
rent carried by the states Both of these quantities 
are affected by the extent of admixture of spin-up and 
spin-down components within the energy eigenstates. An 
initial analysis of some of the results we will discuss has 
appeared previously 

The measurements we describe were performed us- 
ing tunneling devices containing an Al particle less than 
10 nm in diameter, connected to Al electrodes via alu- 
minum oxide tunnel junctions. A device schematic is 
shown inset to Fig. 1(a). The fabrication steps have been 
described in detail previously |^ . An aluminum electrode 
is first deposited on one side of an insulating silicon ni- 
tride membrane containing a 10-nm-scale through-hole. 
The Al is oxidized to form a nm-scale tunnel junction 
near the base of the hole. A layer of Al nanoparticles is 
then formed on the other side of the membrane by de- 
positing 2.5 nm of Al, which balls up into small particles 
due to surface tension. In some of the devices described 
in this paper, the Al evaporation for the particles was in- 
terrupted half-way through and a thin layer of gold was 
deposited to give roughly a 4% (atomic) dose of An inside 
the nanoparticle. Since Al and Au are sufficiently mis- 
cible to form several intermetallic compounds |l4|| , and 
both have significant surface mobilities on the nm length 
scale, we expect that the two types of atoms will be in- 
termixed. When the nanoparticle deposition is complete, 
their surfaces are oxidized to form tunnel junctions, and 
a thick aluminum film is deposited as a second electrode. 
Devices in which tunneling occurs via a single nanopar- 
ticle joining the two leads are selected based on the mea- 
surement of a "Coulomb-staircase" current- voltage curve 
(Fig. 1(a)). 

II. Effects of Spin-Orbit Interactions on Discrete States 

Tunneling spectra of the discrete energy levels are 
shown in Fig. 1(b), at different values of the applied mag- 
netic field, for an Al particle in which we will identify the 
presence of spin-orbit scattering. This particle is nomi- 
nally pure Al, but we have also observed all the features 
that we will ascribe to S-O scattering in Au-doped parti- 
cles. We speculate that the source of the S-0 scattering 
in the nominally pure Al particle is an unintended de- 
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FIG. 1. (a) Large scale Coulomb-staircase curve for a tun- 
neling device containing a nm-scale Al particle at T = 50 mK. 
Inset: Cross-sectional device schematic, (b) Tunneling spec- 
trum of discrete state resonances in the same sample, for a 
range of applied magnetic fields, at T = 50 mK. The curves are 
offset in dl /dV for visibility. Orbital state #2 gives small but 
visible resonances at low B. Small changes in offset charge 
occurred between the 0.1 and 1 Tesla scans and between the 

6 and 7 Tesla scans, shifting peak positions. The 0.1 and 

7 Tesla scans have therefore been shifted along the voltage 
axis, to give the best fit to a linear dependence for peak 1 |. 
The lines tracing under peaks are guides to the eye. 

feet or impurity. Each peak in the dl/dV vs. voltage 
spectrum corresponds to the threshold for electron tun- 
neling via a different quantum-mechanical state in the 
particle, each with the same number of electrons (either 
one more or one less than in the = ground state of 
the particle). In order to convert from the voltage scale 
to the true energy within the nanoparticle, it is neces- 
sary to determine the ratio of the capacitances of the 
particle to the two electrodes, C1/C2. This is measured 
most accurately either by comparing the positions of the 
tunneling peaks due to the same state at opposite signs 
of bias voltage, or by measuring shifts in peak positions 
for superconducting vs. normal-state electrodes The 
conversion factor from voltage to energy for the data of 
Fig. 1(b) is eCi/(Ci + C2) = e(0.53 ± 0.01). A rough 
estimate of the volume of the particle can be made based 
on the capacitances of the particle, determined from the 
spacing between steps in the Coulomb-staircase curve, 
AV = e/Csmaiier = 78 mV. Together with the known 
capacitance per unit area of oxidized aluminum tunnel 
junctions, ^ 50 fF/iim^ p^ , and assuming a particle 
shape that is approximately a hemisphere, we estimate a 
particle radius of approximately r = 3 nm for this device. 
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The peaks in Fig. 1(b) have many features quahta- 
tively similar to previous studies of tunneling resonances 
in pure Al. As the applied magnetic field {B), apphed 
parallel to the plane of the Si^N^ membrane in the de- 
vice, is increased from low-field values, each peak splits 
in two, and the energy difference between these pairs in- 
creases linearly with B at low B (Fig. 2(a)) |^6|. This can 
be understood as Zeeman sphtting of the energies of the 
predominantly spin-up and spin-down states associated 
with each orbital eigenstate. The observation of tunnel- 
ing via both of the Zeeman-split states for the lowest- 
energy tunneling state (#1) indicates that the tunneling 
transition corresponds to a change from an even number 
to an odd number of electrons within the nanoparticle 
Within the uncertainties of the measurement, the 
splitting is symmetric around the low-field resonance en- 
ergy, with little shift up or down for the average of the 
Zeeman-split peaks . This indicates that the effect of 
B on the orbital component of the electronic energy is 
much weaker than on the spin component. This is not 
surprising, due to the particle's small size and disorder. 
Because any real nanoparticle will not have a spherical 
shape or a smooth surface, the orbital angular momen- 
tum of the eigenstates will be quenched to zero in the 
absence of an applied field. In this situation, the cor- 
relation scale which describes the effect of the magnetic 
field on the energy eigenstates is expected to be on the 
order of ^Qy/S/Erh/r'^ , where $0 is the flux quantum, 
d is the mean level spacing, and Eth ~ hvp/i'ir) is the 
Thouless energy scale for a ballistic sample For a 

particle with radius 3 nm the expected field correlation 
scale is approximately 30 Tesla. Since this is much larger 
than the fields of interest in our experiment, throughout 
the paper we will assume that the effect of B on the or- 
bital eigenstates within the particle is negligible, so as to 
concentrate on spin effects. 

There are at least 3 features of the data in Figs. 1(b) 
and 2(a) that differ from typical Al particles. Firstly, 
let us define an effective g factor such that the energy 
splitting between Zeeman-split states is AE — gcSfJ-sB 
(to linear order in B), where /is is the Bohr magneton. 
In over 80% of the nominally pure Al samples we have 
examined previously, (/cff = 2±0.05, which is as expected, 
because S-0 scattering is negligible in pure Al, and the 
free electron (7-factor should apply [|o| . In the sample in 
question, however, t/eff is significantly less, and it varies 
from peak to peak: g^s = 1-84 ± 0.03, 1.68 ± 0.08, and 
1.76 ± 0.05 for the 3 resonances in Fig. 2(a). The second 
difference between this sample and past measurements 
concerns level crossings. In pure Al particles with g- 
factors approximately equal to 2 we have not observed 
departures from linear Zeeman splittings when spin-up 
and spin-down levels corresponding to different orbital 
states cross as a function of B. For a sample without 
S-0 scattering, this must be the case, for then there 
is no coupling between spin-up and spin-down states in 
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FIG. 2. (a) Energies of the discrete electronic states 
within the nanoparticle of Fig. 1, calculated by multiplying 
the voltage positions of the resonances by the capacitance ra- 
tio eCi/(Ci -I- C2) = e 0.53. The thin lines are extensions of 
the low-field linear dependence of the energies on B. Heavy 
lines show the result of the spin-orbit interaction model, de- 
scribing the avoided crossing between levels 1 | and 2 |. (b) 
Magnitude of the current increments contributed by each of 
the first three resonances for positive voltage bias. (Equal to 
the area under the peaks in Fig. 1(b).) Note the crossover in 
magnitude for the current increments associated with states 
1 T and 2 i. 

the Hamiltonian. In contrast, the Zeeman splittings of 
the first two orbital states shown in Figs. 1(b) and 2(a) 
show a clear departure from linear behavior, because the 
upward-trending level from the first state (1 t) under- 
goes an avoided crossing with the downward-trending 
level from the second state (2 |). The third difference is 
that the amplitudes of the resonances for the sample in 
question show unusual features. Whereas in most Al par- 
ticles the amplitudes of the resonances do not display any 
significant S-dependence, here the amplitudes of the two 
levels undergoing the avoided crossing change dramati- 
cally, with the higher-amplitude resonance shrinking and 
the smaller-amplitude resonance growing in the avoided 
crossing region. The amount of current contributed by 
each resonance is plotted in Fig. 2(b), where it can be 
seen that the sum of the current increments contributed 
by these two resonances is approximately constant. 

All of these features can be understood by considering 
the Hamiltonian of the electrons in the presence of S-0 
scattering. The theory behind the physics of the g-factor 
has been considered previously |l^Jl^. Let us write the 
Hamiltonian in zero magnetic field as 



so, 



(1) 
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where Hso contains the terms that couple spin-up states 
to spin-down, and Hq describes all the spin-independent 
forces and interactions. We will neglect the effect of the 
magnetic field on electron orbits, and assume that the 
sample contains no magnetic impurities. Let \An'\) and 
\An\) represent the unperturbed eigenstates of Hq. Then, 
performing perturbation theory to lowest order in Hso, 
the (not normalized) eigenstates of the full Hamiltonian 
have the form 



{Ar^i\fIso\An^)\A^i) 
En ~ Em 



(2) 



The spin-orbit interaction causes the eigenstates to con- 
sist of a linear superposition of spin-up and spin-down 
states; hence the notation ""f" reflects that the eigen- 
state can be considered at most predominantly spin-up. 
The effective g-factor for state n may be written [|lO| 



ffeff,n = ^ 



($„"T"I*«"T") 
\{A^l\Hso\An^)? 



{En — E„ 



(3) 



(Evaluating this expression for the "|" state gives the 
same answer.) The meaning of Eq. (3) is that ges is re- 
duced below the free-electron value of 2 by an amount 
determined by the extent to which S-0 matrix elements 
couple the state n to other states m of opposite spin. Be- 
cause the energy eigenstates are no longer purely spin-up 
or spin-down in the presence of spin-orbit interactions, 
they respond more weakly to an applied magnetic field 
than pure-spin states. Next consider the nature of the 
matrix elements \{Arai\Hso\An^)\'^ . Due to the chaotic 
and strongly fluctuating character of the wavefunctions 
in a metallic nanoparticle Q , the magnitudes of these fac- 
tors will be strongly varying for different values of m and 
n, depending on the details of the wavefunction overlaps 
at the positions of the S-0 scattering defects. Therefore, 
from Eq. (3), it can be seen that different energy lev- 
els in the same sample may have different values of gcs, 
as we observe. Because of the form of the denominator 
in the second term of Eq. (3), we can also expect that 
matrix elements which couple eigenstates nearby in en- 
ergy will produce the strongest influence on geff • We will 
demonstrate an example of this below. 

Let us now begin to analyze the variations in the level 
energies and the currents carried by the particular levels 
displayed in Fig. 2. To do this we will write explicitly 
the form of the effective Hamiltonian matrix for just the 
4 energy levels associated with the first two orbital states, 
which we label as \a j), \a f), \b i), and \b t). The most 
convenient set of basis states arc those which diagonalize 
the spin- independent part of the Hamiltonian, Hq, to- 
gether with all of Hso except that term which couples 
states \a) and \b) to each other. (With this choice, the ba- 
sis states are already not purely spin-up or spin-down, so 



the arrows should henceforth be understood to mean pre- 
dominantly spin-up or predominantly spin-down.) The 
S-0 interaction is invariant upon time reversal. The most 
general Hamiltonian satisfying this symmetry, including 
both ordinary potential scattering and S-0 scattering, 
and describing two Kramers' doublets in the absence of 
an applied magnetic field is (with the above basis choice) 
represented by the matrix : 



H 



(Eai 



d* 
c* 







d 



d 


c 


-c* 


d* 













(4) 



The placement of the zero elements and the arrangement 
of the elements involving c and d are required so that the 
Kramers' doublets are in fact degenerate at B = 0. The 
matrix element d~ {a i\Hso\b i) couples states of the 
same spin, so that it is equivalent to ordinary potential 
scattering for our purposes. Without loss of generality, 
we can pick the orbital basis states \a) and \b) so that 
d = 0. We identify c= {a i\Hso\b T)- Because we are 
assuming that the orbital states are not modified by a 
magnetic field, we take the matrix element c to be inde- 
pendent of B. The only i3-dependence then left in the 
problem is due to the influence of the Zeeman energies in 
the diagonal terms of the Hamiltonian. We write these 
Zeeman energies by including effective g-factors, g'^g, for 
the spin and, simply for convenience in the fitting, we 
also allow a linear term gorblJ^sB (where /u_b is the Bohr 
magneton) to model any shift in the average energy of 
the Zeeman-split pairs. (We will see that the fits give 
(7orti~0.) With these assumptions, the diagonal terms as 
a function of B are 



Ea\.,i — Ea + [gorbM ± 



)fJ'B 



B 



Eb-] ,i — Eb + {gorb,b ± "^^'^ )^-bB. 



(5) 



The terms g'^^^ ^ and g'^g must take into account the S-0 
coupling of state \a) or \b) to all states except each other, 
so that these terms will not be equal to 2. Instead, from 
Eq. (3), we should expect that g^ff ^ and ^^ff b ^^^^ 
related to the total effective g factor by the relationship 



ffeff 



5cff 



\{ai\Hso\b]W 
^ [Ea - Ebf 



(6) 



With d = 0, Eq. (4) gives a very simple Hamiltonian, con- 
sisting of two separate 2-by-2 matrices coupling \a [) to 
l^t), and |at> to |6i>. 

For a weak S-0 interaction, \{a l\Hso\b T)l ^ Eb — Ea, 
the effects of the interaction are easy to understand. 
Away from any degeneracies among the diagonal terms, 
the energy eigenvalues will be approximately equal to the 
diagonal terms, except for a shift in the effective g-factor. 
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When the Zeeman energies arc such that two diagonal 
energies approach degeneracy, they wiU exhibit a simple 
avoided crossing of magnitude equal to 2|(a l\Hso\b T)|) 
because this term couples the two states. Solving the 
Hamiltonian explicitly (with d = 0), the model produces 
an excellent fit for the B-dependence of the measured 
levels (Fig. 2(a)). with the parameters \{a i\Hso\b T)| = 
73 ± 4 neV, g'^s,. = 1-90 ± 0.04, g'^s,b = 1-74 ± 0.04, 
5r„^f, „ = -0.03 ±6.04, and gorb,b = -0.10 ±0.06. The dif- 
ference between the directly measured g- values .goff.i = 
1.84 ± 0.03, 5eff,2 = 1-68 ± 0.08 on the one hand and the 
fitting terms g'^g ^, g'^^ on the other is consistent with 
Eq. (6), since 4| (a j \Hso\b T)lV(^a-Sfc)' = 0.06. From 
this we can see that the S-0 coupling between states \a) 
and \b) contributes approximately 40% of the reduction 
from gcff = 2 for the orbital state 1, and 20% for state 
2. S-0 coupling to other states must account for the re- 
mainder. The fact that we do not have the sensitivity to 
resolve any avoided crossing between the states 2 | and 
3 i (Fig. 2(a)) indicates that the S-0 matrix element 
coupling these states is smaller than (a i\Hso\b t)- 

The changes in the amount of current carried by the 
resonances (Fig. 2(b)) can be understood by examining 
the manner in which the energy eigenstates are composed 
of linear superpositions of basis states. Consider the two 
energy eigenstates {\lower) and \upper)) formed from su- 
perpositions of the avoided-crossing basis states \a t) and 
1 6 I). By diagonalizing the Hamiltonian (Eq. (4) with 
d = 0), it is simple to demonstrate that these have the 
form 



\lower) = -f{B)\a ^) + T]{B)\b I) 
\upper) = -r]*{B)\a^)+riB)\bl) 



(7) 



where the coefficients jiB) and ri{B) depend on B as 
shown in Fig. 3(a). The key point is that, as the mag- 
netic field is varied in the avoided crossing region, the 
relative contributions of \a t) and \b J,) to each eigenstate 
will change, and consequently the tunneling currents can 
be altered. This simple conclusion will be the topic of the 
next two pages of discussion. The reason for an extended 
analysis is that the magnitudes of the currents are de- 
termined by a process of sequential tunneling across the 
two tunnel junctions in the device, so that the measured 
current values are not simply a function of a the tunnel- 
ing rate into an individual energy eigenstate. Instead, 
the current will be affected by all energetically-allowed 
transitions within the device. In order to deal with this 
complication, the plan of our discussion is that we will 
focus first on the bare tunneling rates TL,iower, ^L,upper, 
^R.iower, and Tji^upper for tunneling of an electron be- 
tween the energy eigenstates {\lower) and \upper)) and 
the left (L) and right (R) electrodes. Later we will exam- 
ine two different limits for calculating the total current 
through the device in terms of these bare tunneling rates. 
In either case we will see that, despite the complications, 
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FIG. 3. (a) Magnetic-field dependence of the coefficients 
in Eq. (9), for the superposition of predominantly spin- up 
and spin-down basis states occurring in the avoided crossing 
of levels 1 t and 2 J,, (b) and (c) Markers: Tunneling rates 
for the energy eigenstates in the avoided crossing region, esti- 
mated as discussed in the text using either Eq. (11) or (12), as- 
suming that the relaxation rate of non-equilibrium excitations 
within the particle is either slower or faster than the tunnel- 
ing rates. Lines: Fits using the predictions of the spin-orbit 
Hamiltonian. Regardless of the energy-relaxation rate, the 
magnetic-field dependent changes in the currents flowing via 
the levels in the avoided crossing region (Fig. 2(b)) can be 
explained qualitatively by the change in tunneling rates ex- 
pected from the S-0 Hamiltonian. 



the changes in the total current carried by a particular 
tunneling resonance as a function of magnetic field can be 
related to the changing composition of the energy eigen- 
states in the avoided-crossing region (Eq. (7)). 

For the high-resistance barriers used in the experiment, 
the bare tunneling rates between either of the electrodes 
and energy levels in the nanoparticle can be written in 
terms of matrix elements of a tunneling Hamiltonian Ht 
which couples states in the electrodes to the energy eigen- 
states. Since tunneling of a spin-up electron from the 
electrode is necessarily incoherent with respect to tunnel- 
ing of a spin-down electron, we have for the left junction 
(for the right junction the equations are similar): 
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T X! \\{'^electrodeA\HT\loWer)\'^ 

ie 

+ 1 {ipelectrode.l | i?T | ^OWCr ) ^ | 



LAower 



electr 



(8) 

TL,upper = ^ ^ \^\{i'electrode,]\HT\upper)\'^ 

elecirode 



+ 1 {tpelectrode,l \ Ht \ Upper) ^ | . 



These expressions can be given in a more illuminating 
form by writing llower) and \upper) explicitly as linear 
superpositions of the basis states \a 1) and \b |) (as in 
Eq. (7)). The tunneling rates become 



(9) 



where the -B-independent tunnel-coupling strengths for 
the basis states are 



= ^E4,\{^electrodeMHT\a^)\^ 
— \{^electrode,l\HT\b 



(10) 



For values of B well below the avoided crossing range, we 
have j{B) « 1, r]{B) 0, and the rates for tunneling into 
the energy eigenstates are equal to T^^ai and f,x- Since 
in this regime the total current passing through the \b) 
resonances (orbital state #2 in Figs. 1(b), 2(b)) is very 
small compared to the \a) peaks (orbital sate #1), clearly 
these two rates must be very different. As B is swept 
through the avoided crossing region, the admixture of 
the two basis states within the eigenstates changes, with 
17(5)1 evolving gradually from 1 to 0, and |?7(-B)| going 
from to 1. This means that there should be a grad- 
ual exchange of tunneling weight between \lower) and 
\upper), with TL.iower evolving from T^.a] to TL,bu and 
^L,upper doing the reverse. Qualitatively, this crossover 
behavior is apparent in the currents in Fig. 2(b). 

In order to attempt a more quantitative treatment of 
the measured currents, it is necessary to analyze the re- 
lationship between the bare tunneling rates F discussed 
above, and the value of the current that results from 
sequential tunneling across the two tunnel barriers. This 
requires a full consideration of all the processes that 
can occur during current flow. When the applied bias 
is larger than the level spacing, non-equilibrium elec- 
tron distributions are produced within the nanoparticle 
during tunneling, and these can open new channels for 
electrons to flow []6|,pO|. The idea is shown in Fig. 4. 
In Fig. 4(a) we show the simple process of an electron 
tunneling from the left electrode to an empty level on the 
particle, when the voltage across the device is sufficient 
to supply the threshold tunneling energy. Due to electro- 
static interaction with this additional electron upon its 
arrival, the lower-energy electronic states already filled 
within the particle can be shifted up in energy to the 
positions drawn. If the applied voltage needed to initi- 
ate tunneling is larger than the level separation between 



(a) 


— •— 




(b) 


— • — 
• 










— •— 






— • — 




> 






— •— 













(c) 












— • — 






FIG. 4. (a)-(d) Some of the allowed transitions contribut- 
ing to the magnitude of the current flowing at a tunneling 
threshold. 

states, one possibility for the next step in the tunneling 
process might be as shown in Fig. 4(b), where an elec- 
tron tunnels out of one of these lower-energy states to 
the right electrode, leaving an electron-hole excitation on 
the particle. After this, the excited electron might relax 
(Fig. 4(c)) before the next electron tunnels onto the par- 
ticle. Alternatively an electron may tunnel from the left 
electrode into the hole (Fig. 4(d)). All of these processes, 
and all other energetically- allowed transitions, will con- 
tribute to the value of the current that results when the 
voltage is turned above the threshold required to initiate 
tunneling in Fig. 4(a). The lower-energy, initially-filled 
states can contribute to the current even though they are 
hidden in the sense that they do not produce tunneling 
thresholds of their own. In general, to calculate the to- 
tal current, one must solve a full master equation which 
takes into account the rates of all the allowed transitions, 
including the hidden levels. The parameters entering the 
calculation are the number of hidden levels, the tunneling 
rates between each level and the left and right leads, and 
the relaxation rates for the different allowed excitations 
within the particle. 

For samples of the sort investigated in Figs. 1-3, which 
do not have a gate electrode that can be used to ad- 
just the number of hidden states, the measured data are 
not sufficient to fully determine all of the parameters re- 
quired to describe the currents quantitatively. However, 
progress can be made with some simplifying assumptions. 
One particular difficulty is that we do not know how the 
relaxation rate for excitations within the particle {e.g., 
Fig. 4(c)) compares to the tunneling rates in this sam- 
ple. The predicted order of magnitude for the relaxation 
rate due to electron-phonon scattering is lO^s^^ ||^, less 
than the tunneling rates we will determine, but only by 
about a factor of 10. Allowing for some uncertainty in 
the theory, we will consider both of the two simple limits 
- that the relaxation rate is cither much slower or much 
faster than the tunneling rate. The following discussion 
should not be considered a quantitative determination 
of tunneling parameters, but it will serve to illustrate 
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the way in which the measured changes in current in- 
crements as a function of B can be hnked to the bare 
tunneling rates. We also invoke three other simplifying 
assumptions: (1) that the number of hidden states stays 
the same over the range of voltage and magnetic field an- 
alyzed in Figs. 1(b) and 2, (2) that the ratios Ti^i/Tji^i, 
for the tunneling rates from quantum level i to the left 
and right electrodes, are all the same, and (3) that for 
all hidden states the tunneling rates to the left electrode 
are identical = r^^h- We define x = TL/Tii. These as- 
sumptions reduce the free parameters in the problem to 
a tractable number, but we make no rigorous claims as 
to plausibility. 

We first consider the slow-relaxation limit, in which we 
can ignore all processes of the sort pictured in Fig. 4(c). 
In this limit, the second assumption listed above leads to 
a great simplification, because the probability for occupa- 
tion of any quantum level accessible by tunneling will be 
the same. If the voltage bias is such that N hidden levels 
and M originally-empty levels participate in tunneling, 
the total current that results at T = is, by solution of 
an elementary set of rate equations, for positive (+) and 
negative (— ) bias: 



Im. 



N,+ 



lM,N,+ 



lM,N,- 



eM{N + l){NTL,h + Y.Zl^L,^) 
{N + M){N + I + xM) 

eM{N +l){NTL,h + Y.ti^L,i) 
x{N + M){N + 1 + 



(11) 



Within the model, we can determine the values for the 
4 free parameters N, a;, F^ /j, and F^ i (the tunneling 
rate into the lowest-energy initially-empty orbital state 
at low field) from the four low-field {B — 2 Tesla) values 
of the total, cumulative currents flowing at the 1 t and 
1 I thresholds, for positive and negative bias: /i,jv,+ = 
A/i4,+ = 9.6 X 10-" A, I2.N.+ = A/1T.+ + A/1X.+ = 
1.41 X 10-1" A, h,N,- = A/ix,_ = 6.6 X 10-" A, 
and h^N- = A/11 _ + A/ix,- = 1.19 x 10-^° A. Be- 
cause of time reversal symmetry we can assume that 
Tl.iT — ^L,ii for small B. The results for the four free 
parameters are N = 2.4, x = 2.0, F^^/j = 9.4 x 10*^5-^ 
and Fl_i = 1.01 x 10^s~^ The fact that N is not an in- 
teger may reflect the weaknesses of the assumption that 
the ratio x = T l /T r is the same for all the quantum lev- 
els and/or the assumption of slow relaxation. Employ- 
ing these values and the measured (positive-bias) cur- 
rent increments for the avoided-crossing states shown in 
Fig. 2(b), we can then invert Eq. (11) (for M =2 and 3) to 
estimate the bare tunneling rates TL^io^er and F^. upper- 
over the whole range of B from 1 to 7 Tesla, with the 
results shown in Fig. 3(b). 

In the same way we can consider the fast-relaxation 
limit, in which the electrons in the nanoparticle relax to 
their lowest energy state between all tunneling events. 
The solutions to the rate equations are 



(12) 



In this case there are just three parameters, F^ i, x, and 
Nr^^h, which can be determined from the B — 2 Tesla 
values of /i,Ar.+ = A/i|^+, l2,N,+ — A/i|^+ + A/i|_+, 
and h^N- = ^hi - as Fl,i = 9.3 x 10^s-\ x = 1.96, 
and NTlji — 2.4 x lO^s-^. These parameters, together 
with Eq. (12) predict a value of 1.18 x 10-^" A for 
I2.N,- — A/i|^_ + A/i|__, in good agreement with the 
measured value, 1.19 x 10"^" A. We can then invert 
Eq. (12) using the measured current increments of the 
avoided-crossing states in Fig. 2(b) to estimate the tun- 



neling rates 



lower 



and Fl 



upper 



in the fast relaxation 



limit (Fig. 3(c)). The differences between Figs. 3(b) and 
3(c) reflect to some extent the degree of uncertainty with 
which we can estimate these bare tunneling rates. 

We see from both Figs. 3(b) and (c) that the crossover 
observed in the magnitude of the current increments 
for the two avoided crossing states can be related to a 
crossover in the bare tunneling rates, of the type pre- 
dicted by the spin-orbit scattering Hamiltonian. The 
lines in Fig. 3 (b,c) display fits to the S-0 model result, 
Eq. (9), with only 2 adjustable parameters, T^^a'i, and 
TL.bi, which simply set the B = values of the tunneling 
rates. For the slow-relaxation limit the B-dependence 
of the tunneling rates is very well described by the S- 
O formalism. In particular, as predicted by the model, 
the tunneling rates cross close to the same magnetic field 
value, 5.4 Tesla, where the avoided-crossing levels have 
their closest approach. Also, the tunneling rates well be- 
yond the crossover regime are approximately equal to the 
B = tunneling rates. Neither result holds for the cur- 
rent increments themselves (Fig. 2(b)), due to the effect 
of the hidden levels. The agreement between the S-0 
theory and the tunneling parameters estimated in the 
fast-relaxation limit is not quite as close as for the slow- 
relaxation limit. This is consistent with the estimates in 
ref. that the energy relaxation rate is 10^5"^, an order 
of magnitude less than the tunneling rates we determine. 



III. Effects of Spin-Orbit Interactions on 
Superconducting Nanoparticles 

We next consider different samples, a larger Al particle 
containing 4% Au impurities (Fig. 5(a)), compared to a 
pure Al particle of similar size showing no indications of 
S-0 scattering (Fig. 5(b)) Because of their larger 

size, the mean level spacings in both samples are smaller 
than in the device of Figs. 1-3, but nevertheless a large 
energy difference is visible between the first and second 
peaks in both spectra. This is characteristic of odd-to- 
even tunneling in a superconducting particle. The energy 
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FIG. 5. (a) Tunneling spectrum for an Al particle con- 
taining 4% Au impurities, for a sequence of magnetic fields 
from 0.03 to 9 Tesla in 1 Tesla increments, T = 15 mK. The 
particle exhibits a superconducting gap for odd-to-even elec- 
tron tunneling, (b) Comparison data for a pure Al particle, 
from ref. [5]. The curves are artificially offset for visibility. 

gap reflects the large difference in energy (approximately 
twice the superconducting gap A) required for the tun- 
neling of an electron to reach the ground state of an even- 
electron superconductor in which all electrons are paired, 
versus the first excited state containing two unpaired 
quasiparticles Q . The two samples display qualitatively 
different behavior in a number of respects, however. One 
difference is that the resonance peaks in the Au-doped 
sample are somewhat broader. We believe that this is 
not related to the impurities, but is instead an effect of 
non-equilibrium distributions of electrons on the island 
1^,0, excited by a source-drain voltage 3 times as larger 
than what is needed to overcome the Coulomb blockade 
and initiate tunneling in Fig. 5(a). Similar broadening 
could also be observed for the particle in Fig. 5(b) (see 
ref. @, Fig. 4), when a gate voltage was used to shift the 
tunneling spectrum to comparable values of the source- 
drain voltage. The device in Fig. 5(a) had no gate. 

We will focus instead on the differences in the magnetic 
field dependence of the data in Fig. 5(a) and (b). For 
the pure Al superconducting nanoparticle, the primary 
effect of a magnetic field is to produce linear shifts corre- 
sponding to a Zeeman-spin-splitting with g^g = 2 ± 0.05 
p2[ . The gap in the spectrum decreases linearly due 
to the difference in Zeeman energies for the ground and 
first-excited states of the superconductor, until it goes 
to zero at about 4 Tesla. Models of superconductivity 



in small particles |^,^ relate this crossing with the su- 
perconducting critical field, because the extent of elec- 
tron pairing correlations drops abruptly at this point, 
although fluctuation-induced effects of attractive elec- 
tron interactions may persist In contrast, all the 
resonance energies for the sample containing gold impu- 
rities are significantly less sensitive to an applied mag- 
netic field. Instead of Zeeman splitting with ^eff = 2, 
the ground and first-excited-state transitions move at 
low B with slopes 5eff,i/2 -I- gorb,i = 0.41 ± 0.03 and 
— 5cff,2/2 + gorb,2 = —0.27 ± 0.03, suggesting values for 
geS in the range 0.5-0.8. Even at 9 Tesla, the gap between 
these states has not gone to zero, indicating a much larger 
critical field for superconductivity in this sample than for 
pure Al. Similar increases in critical fields due to the re- 
duction in the effective g-factor caused by S-0 scattering 
are familiar for thin films in parallel magnetic fields, and 
in other contexts where superconductivity is limited by 
spin-induced pair-breaking |26|. At fields above 6 Tesla, 
the slope of the energy vs. B curve of the ground-state 
transition in the Au-doped sample changes sign (with the 
energy decreasing with increasing B at high fields), sug- 
gesting an avoided crossing with the higher- lying levels. 
The minimum gap between the ground-state and first 
excited-state peaks corresponds to a S-0 matrix element 
of magnitude approximately 130 fieV p7|] . 

The presence of S-0 scattering must necessarily change 
the nature of the superconducting transition in a mag- 
netic field. As we noted above for pure Al particles, 
the extent of superconducting pairing correlations is pre- 
dicted to drop abruptly at the magnetic field for which 
the energy of the first state that moves to lower energy 
with increasing B (meaning that it is a spin-1 tunneling 
state) crosses below the energy of the upward-trending 
(spin-0) ground state, so that it becomes energetically 
favorable to break a Cooper pair. In contrast, in the 
particles with significant S-0 scattering, the existence of 
avoided crossings means that energy levels correspond- 
ing to different spin states do not cross. Therefore the 
disruption of pairing correlations must occur gradually, 
as the spin content of the particle's ground state changes 
continuously in the avoided crossing region. 

Notably, the magnitude of the superconducting gap at 
B = is not significantly affected by the presence of S- 
O scattering. Setting the difference between the ground 
and first-excited state energies equal to 2A, we have A w 
0.25 meV for the pure Al particle of Fig. 5(b), similar 
to previous values and A « 0.26 meV for the par- 
ticle with Au impurities. This similarity is as expected, 
since S-0 scattering does not break time-reversal symme- 
try and therefore does not interfere with superconducting 
pairing. 
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IV. Conclusions 

We have examined a number of effects associated with 
the presence of S-0 scattering in metal nanoparticles. 
The sensitivity of the "electron-in-a-box" energy levels 
to an applied magnetic field is weakened, so that they 
exhibit effective t;- values less than 2. When predomi- 
nantly spin-up and spin-down levels approach each other 
as a function of magnetic field, they may undergo avoided 
crossings due to S-O-induced coupling between the spin- 
up and spin-down states. In the avoided crossing region, 
the magnitude of the current transported at the reso- 
nance thresholds can change, due to the changing ad- 
mixture of spin-up and spin-down basis states that com- 
prise the energy eigenstates. The presence of Au im- 
purities does not greatly modify the size of the super- 
conducting gap in Al particles large enough to exhibit 
superconductivity. However, the critical magnetic field 
for the destruction of superconductivity is increased. We 
also argue that with the presence of S-0 scattering, the 
superconducting pairing parameter should vary contin- 
uously at large fields, because S-0 scattering eliminates 
the simple level crossings which cause the extent of pair- 
ing correlations to drop abruptly in pure Al samples. As 
a final remark, we note that all of the results that we de- 
scribe can be adequately explained by treating the S-0 
interaction perturbatively, and by ignoring the effect of 
the magnetic field on electron orbits. For samples with 
stronger S-0 interactions or with larger sizes such that 
effects of an applied field on the orbital states are signifi- 
cant, a more sophisticated treatment would be necessary 
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